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Abstract 

Polynomial Lie (super)algebras g p d are introduced via Gj-invariant polynomial 
Jordan maps in quantum composite models with Hamiltonians H having invari- 
ance groups Gj. Algebras g p d have polynomial structure functions in commutation 
relations, are related to pseudogroup structures expV, V £ g p d and describe dy- 
namic symmetry of models under study. Physical applications of algebras g p d in 
quantum optics and in composite field theories are briefly discussed. 

1. INTRODUCTION 

The symmetry methods are fruitfully exploiting in quantum physics from the time 
of its origin (see, e.g., [1-6] and references therein). In particular, they provide an 
elegant treatment of physical tasks using a powerful formalism of Lie groups and alge- 
bras, especially generalized coherent states (GCS) and related technigues [2,4-6], when 
Hamiltonians H of systems under study are linear forms 

H = KF a + C, [F at F P \. = F a F/3 - FpFa = £ c^F 7 ee ^({F 7 }) (1.1) 

a=l 7=1 

in generators F a of d(< oo)-dimensional Lie algebras g D of dynamic symmetry (A a are 
c-number coefficients and [F a , C] = 0). But for last years nonlinear models in many 
branches of quantum physics have called for different extensions of usual Lie algebras 
in Eq. (1.1) by means of 1) admitting infinite dimensions d [7,8], 2) involving two types 
commutators ([, ] T ) defining Lie superalgebras [7] and 3) using nonlinear structure func- 
tions ip a/ 3 {{F a }) defining nonlinear or deformed Lie (super) algebras [7-11]. However, 
at present, these new Lie-algebraic structures, enabling to display important structure 
features of models with such generalized Hamiltonians (1.1), do not yield universal tech- 
niques for solving physical tasks due to the absence of simple ("finite") "disentangling" 
(Zassenhaus) and "multiplication" (Baker-Campbell- Hausdorff) formulas for their ex- 
ponentials [5,7,10] which determine a high efficiency of group-theoretical methods [2,4]. 
Therefore, keeping in mind relevant extensions of these methods for new algebras, it 
is of importance to examine possibilities of adequate modifications of group-theoretical 
and Lie-algebraic techniques for single classes of models. 



In the present paper, summarizing and developing results of the papers [5, 12,13], 
we discuss these problems for polynomial Lie (super) algebras g p d which describe dy- 
namic symmetries of composite many-body models with Hamiltonians H having invari- 
ance groups Gi(H) and are obtained via generalized Gj-invariant Jordan maps [5]. 

2 . POLYNOMIAL (^-INVARIANT JORDAN MAPS AND POLYNOMIAL 
LIE (SUPER) ALGEBRAS IN MANY-BODY PHYSICS 

As is known, in composite models of many-body physics whose Hamiltonians 
H and quantum state spaces L(H) are given in terms of boson-fermion operators 
cii, af, bj, bj~([a,i, alj}^ = 6ij,[bj,bj~] + = Sij), Lie-algebraic methods are introduced in 
a natural way via different boson-fermion maps [2-5,14]. Specifically, the Jordan map 
[5,14]: (di,af,bj,bj~) i— > F a giving generators F a by quadratic forms in a,i,af,bj,bj 
reduces quadratic (in field operators) Hamiltonians H (ai,af ,bj,bj) to the form (1.1) 
[2-5]. 

This map, introducing collective dynamic variables F a e gff, is particularly fruitful 
when Hq have (both continious and discrete) invariance groups Gi(H ): 

[Ho, Gi(H )]_ = [£ D , Gi(H )]_ = (2.1) 

and field operators are transformed with respect to certain (as a rule, fundamental) irre- 
ducible representations (IRs) of the groups Gi(H ) [5,15]. Then F a e g$ are quadratic 
vector Gi(Ho)- invariants; besides, by virtue of the construction and Eq. (2.1), invari- 
ant (Casimir and class) operators Ck{Gi) and Ck(go) determining IRs of Gi and g® are 
functionally connected and their eigenvalues on spaces L(H) are specified by certain 
common sets [U] = [l , li, . . .] of invariant quantum numbers labeling extremal (lowest) 
vectors \ of both Gi and g® IRs. All that, in turn, entails spectral decompositions 

°([li])W]) (2-2) 

Ik] 

of spaces L(H) in direct sums of the subspaces L ([li]) which are invariant with respect 
to joint actions of algebras gff and groups Gi being carrier-spaces of so-called factor- 
representations (isotypic components) [2] of both algebraic structures Gi and gff. 

In the case of suitable groups Gi decompositions (2.2) have the simple spectra 
(cr([/j]) = 1), and, then, pairs (G^g®) say to act complementarily [15] on L(H) and 
to form the Weyl-Howe dual pairs since such pairs were first considered in quantum 
mechanics by H. Weyl in the analysis of interrelations between unitary and permutation 
symmetries of iV-electron systems [1], and their explicit mathematical characterization 
was given by R. Howe [16]. A physical importance of dual pairs is due to that they 
describe completely both invariance and dynamic symmetries of models under study 
(see, e.g., [15,5,12] and references therein). 

The constructions above are generalized in a natural manner when extending 
quadratic Hamiltonians H by (^-invariant polynomials Hj(ai,af ,bj,b~j~) of higher de- 
grees which describe essentially nonlinear interactions [5,17]. Then generalized dual 
pairs (Gi,g p d) of invariance groups Gi and (describing dynamic symmetry) Lie-like (su- 
per) algebras g D = g p d are obtained via (^-invariant generalized Jordan maps [5] 

(a i ,a+,b j ,bf)^(F a ,V x ,V+) (2.3) 

expanding the sets {F a e g^} by some additional generators V\, which are si- 
multaneously elementary vector Gi-invariants and components of two mutually con- 
tragradient g$- irreducible tensor operators V, V + given by homogeneous polynomi- 
als in ai,af,bj,bj. (Note that in practice Hamiltonians H ,Hj may contain, besides 



Oj, af, bj, b~j~, other -covariant operators, e.g., the Pauli matrices cr a (i) etc. that leads 
to appropriate modifications of Eq.(2.3) [5,12].) Then, by virtue of the vector invariant 
theory [1], the sets {F a ,V\} form finite-dimensional integrity bases [1] of associative 
algebras of Gj-invariants embedded in enveloping algebras U(w(m)) of the Weyl- 
Heisenberg (super) algebras w(m) with generators di,af, bj, b~j~. Furthermore, endowing 
sets {F a ,V\} by commutators [, ]±, one gets (via a specific extension of the Ado's 
theorem [2]) finite-dimensional Lie-like (super) algebras g p d = g D of dynamic sym- 
metry which, however, have polynomial structure functions ^^({-^qJXp = deg{ip)) for 
commutators [V a: Vg]±, [V a: Vp]± (the subscripts "±" are determined by the "fermion 
contents" of operators V a ) and may be named as polynomial Lie (super) algebras. Em- 
phasize an importance of the Gj-invariance of polynomials Hj(. . .) because, in general, 
it is impossible to get finite-dimensional algebras if cancelling this condition [5,12]. 
Algebras g p d are extensions of Lie algebras gff and have the coset structure [10,11]: 

9 P d = h + v, h = g°, v = Span{V a , V+}, [h, v] C v, [v, v] C U{h) (2.4) 

that enables us to construct IRs of g p d starting from /i-modules [5] (U(h) are enveloping 
algebras of h, and hereon we omit the subscript "±" in [, ]±). For example, exten- 
sions via (2.3) of the unitary algebras u(m) by their (G„-invariant) symmetric and 
(SU (n)-invariant) skew-symmetric tensor operators give two classes of polynomial os- 
cillator (super) algebras (see Section 4) wheras such extensions of the symplectic algebras 
sp(2m,R) by S'0(n)-invariant skew-symmetric tensors yield polynomial deformations 
of the Lie algebra u(m,m) [5]. Without dwelling on other general properties of al- 
gebras g p d we only note that they are close in their structure with W n -algebras [10], 
can be enlarged (via repeated commutators) to certain graded infinite-dimensional Lie 
(super) algebras g p d = Yl^LoQr, [<7rj<7s] C g r +s, 9r>o = W(/i)(t>) r , and their exponentials 
exp(g pd ) generate (non-analytical) pseudogroup structures having, in general, no finite 
disentangling formulas [5]. And now we consider applications of concrete algebras g p d- 

3 . POLYNOMIAL LIE ALGEBRAS sl pd {2) IN QUANTUM OPTICS 

Simplest examples of polynomial Lie algebras are given by algebras s/ p ^(2) = 
Span(V , V±), obtained via extending the unitary algebra u(l) = Span(V ) by gen- 
erators V± and satisfying the commutation relations (CRs) 

[V , V ± ] = ±V±, [V-, V + ] = r~\V Q ) = V p (V + 1) - ^ p (Vo), [V p (Ro), V a ] = (3.1) 

where \I /P (. . .) is the polynomial of degree p in the variable Vo, ^ p (Rq) = \l/ p (Vo) — V + V- 
is the sl p d(2) Casimir operator (with R being the "lowest weight operator") and hereon 
we omit the identity operator symbol / in expressions like ^ p (y + al). As is seen from 
Eq. (3.1), algebras sl pd (2) are reduced to the Lie algebra si (2) when p = 2 and may be 
considered as its specific deformations; they are also obtained from certain g-deformed 
algebras by means of the Wigner- Indnu contraction when q — > 1 [6]. 

Algebras sl pd (2) arise via the map (2.3) in nonlinear models of quantum optics 
where coset generators V± are interpreted as creation/ destruction operators of specific 
coherent structures (clusters) [12]. For example, Hamiltonians 

H = H + H I = j2 ^afa, + g{a+) n {a ) m + g* (a+) m (a x ) n , m < n, (3.2) 

i=0 

describing multiphoton processes of scattering on multimode Fock spaces Lp = L(H) = 
Span{\{rii}) = ^[^.^^(af^lO >} (g are coupling constants, uji are field modes 



frequencies and H = 1), have invariance groups Gi(H) = C n <8> C rn <S> exp(i\Ri) where 
C n = {ckn = exp(i2nk/n) : af — > Ck n af},Ri = (mafai + na^ a ) / (m + n) . 
Then the map (2.3) given as follows 

V = (a+ ai - a+a )/(m + n), V + = (a+)"(a ) m , V. = (V + ) + (3.3) 

reduces Eq. (3.2) to the form 

H = aV + gV+ + g*V- + C, [V a , C\ = 0, a = nuj l - mu , C = R^ + u ) (3.4) 

and determines the generalized dual pairs (Gi(H), sl p d(2)). The structure polynomials 
ty p (V ) are determined with the help of Eqs. (3.3), the characteristic relation 

(y + V--**(y ))\ L{H) = o (3.5) 

and defining relations for a(i),a + (i) [12]: 

*P(Vb) = {nV +Ri) (n) {Ri-mV +m) {m \ p = m+n, A [b) = A{A-l)...{A-b+l) (3.6) 

(An extra dependence of \& p (Vo) on R\ reflects functional interrelations between in- 
variant operators of Gi(H) and sl p d(2).) The subspaces £([^]) = Span{\[li\; v) =oc 
Vv\[h\), V \[h];v) = (k + v^v),^^}^) = h\[h};v),i = 0,1, V-\[U]) = 0} in 
Eq. (2.2) are generated by the lowest vectors |[Zj] >= [s!/t!] _1 / 2 (ao ) s (a^ ) K |0 >, k = 
0, 1, . . . , n — 1, s = 0, 1, ... where l = (k — s)/{m + n), l\ = (ttik + ns)/(m + n); k 
specifies IRs of discrete invariance subgroups C n and s = d([li\) for compact versions 
of sl p d{2) that is the case when m ^ in (3.2). (Compact (su p d(2)) and non-compact 
(su p d(l, 1)) realzations of sl p d(2) algebras are distinguished depending on whether di- 
mensions d([li\) of the spaces £([/i]) are finite or infinite [12].) Eqs. (3.1), (3.4)-(3.6) 
yield requisites for developing both exact and approximate Lie-algebraic methods to 
solve physical tasks in models (3.2). We outline them following the papers [12,13]. 

Exact methods are based on using Eqs. (3.1) and their resemblances with defining 
relations for si (2). Thus, substituting Eq.(3.4) in the Heisenberg equations for cluster 
dynamic variables V a (t) related to generators V a one gets non-linear equations 



dV T/ _ ,dV+ „ .dV- 



gV + - g*V_, i-± = -aV + - g*iT\V Q ), i— = aV. + g^-\V Q ) (3.7) 



dt dt dt 



generalizing linear Bloch equations for sl(2) and having in the cluster mean-field ap- 
proximation (< I >= /(< \V a \ >) quasiclassical solutions expressed in terms of 
hyperelliptic (Abelian) functions [12]. However, such direct extensions of s/(2)-algebraic 
techniques [4] are impossible for finding evolution (Unit)) and diagonalizing (S) oper- 
ators because relevant disentangling formulas forexp(J2i a-iVi) and explicit expressions 
for matrix elements ([/«]; f\ exp(J2i a>iVj)\ [/»]; v) are absent [12]. 

Nevertheless, substituting " pseudogroup" (cf. [9]) representations 



Uu(t) = £ V{uf(V ;t), S= £ V{S f (V ) (3.8) 

f=—oo f=—oo 

for U H {t) and S (with V~ k = V* [nf=o ^ P ( V o ~ Op due to Eq. (3.5)) or expansions 



i^([y;/)) = ^/n(^- A J(^))iw)=^n(^-4)iw) = E^(^/)iM^) (3-9) 

3 3 v 



for energy eigenstates \E([li\; f)) in the diagonalizing scheme SHS = H(Vq) and the 
time-dependent Schroedinger equation ihdU u{t) / dt = HU(t) and using Eqs. (3.1), 
(3.5), one gets finite-difference and differential- difference equations determining (to- 
gether with unitarity conditions) "coefficients" Sf(Y ), u^(V ;t), diagonal Hamiltonian 
forms H(V ), amplitudes Q v (Ef) and energy spectra {£([/;]; /)} [12,13]. (Note that two 
first equalities in (3.9) realize, in fact, the algebraic Bethe ansatz [17] for wave functions 
\E([li\; f)) in terms of the sl p d(2) generators [13].) These equations define new (non- 
classical) orthogonal functions in both discrete and continious variables which are simul- 
taneously related to solutions of singular differential equations yielded by using two con- 
jugate differential realizations of generators V a [12,13]: V + — z, V — zd/dz + l , V- = 
z~ 1 ^{zd/dz + l ) and V_ = d/dz, V = zd/dz + l , V + = ^(zd/dz + ^(d/dz)- 1 gen- 
eralizing the Bargmann or GCS representations [2,4] for sl(2) generators. However, 
at present, simple analytical expressions for these special functions are absent in gen- 
eral cases though some integral representations were found for them with the help of a 
specific "dressing" of s/(2)-solutions of certain auxiliary exactly solvable tasks [12]. 

Approximate methods developed in [12,13] are based on realizations of generators 
V a as special elements of extended enveloping algebras Uq,{sl{2)) of the si (2) algebra 
via a generalized Holstein-Primakoff map given on each subspace £([/«]) as follows [12] 

V = Y + l ± J, V + = y + [<£(Yo)] 1/2 , *(Y ) = P^Yy V- = (V + )+ (3.10) 

where Y a and ^ 2 (Y ) — (J ± Y )(±J + 1 — Y ) are generators and the structure poly- 
nomials of the su(2) / su(l, 1) algebras, =fJ are lowest weights of their IRs realized on 
L([/j]) (upper/lower signs correspond to su p d(2)/su p d(l, 1) versions of sl p d{2)). 

Eqs. (3.10) enable us to re-write restrictions of Eqs. (3.4) on £([/«]) and basis 
vectors |[Zj]; v) in terms of Y a : 

H [li] ({Y a })=aY + gY +y /^j+g* y /^jY.+C, \[% v) = M(J, v)(Y + y\[k]) (3.11) 

(with C = C([li]) + a(l ± J)) and to use the formalism [4] of the SL{2) GCS 

\U,v;0 = S Y \[%v) = J2Sj v (0\UJ), S Y (0 = exp(er + - CY-) (3.12) 

f>o 

(with Sj v (£ = re 10 ) being the SL(2) Wigner D-function) for examining models (3.4) in 
SX(2)-cluster (taking into account mode correlations in (3.2)) quasiclassical approxi- 
mations [13]. Specifically, SL{2) GCS determine quasiclassical energy functionals 

C + a{v T J)c(2r) - [ge^ + g*e»\ E/>o 1^(0 WO Wo + 1 + f)] 1/2 1 J 
(with c(r) = cosr/ coshr for su{2) / su{l, 1)) or their mean-field approximations 

w cm/ (W;«;0 = ^M«Pi];«;eiTOIW;«;0)) = 



C + a(v =F J)c(2r) - [ge~ ie + g*e ie ](J =f v)s(2r)[$((^J + v)c(2r))\ 



i/ 2 (3-14) 



(with c(r) = cosr/ coshr, s(r) = sinr/ sinhr for su(2)/su(l, 1)) which can be applied 
in standard calculation schemes [2,4,13]. 

For instance, inserting Eqs. (3.13) and (3.14) in the stationarity conditions 

a) gM!Mpi) =0 , t) W(M;.;f) =0 (315) 



one finds approximate eigenfunctions \E cq / cm f ([li]; v )) = SV(£o)^| [h]', v) and appropriate 
eigenenergies E<*([%v) = H c %[h};v^ ), E cm f([h};v) = H cm f ([% v; where £ = 
r o9/\9\ an d values r depend on g, a, U and are determined by real solutions of algebraic 
equations obtained from Eq. (3.15b) with v — [13]. 

The energy functionals (3.13), (3.14) may be also used for a quasiclassical analysis 
of the SX(2)-cluster dynamics described by the classical Hamiltonian equations (cf. [4]) 

q=^, p= ~l£> n={z(t);[l i ]\H\[l i ];z(t)) (3.16) 

for the canonical parameters p = (z(t); [li]\Y \[li]; z(t)) = =FJc(2r),g = 6 of the SL(2) 
GCS |[Zj];z(t) = rexp(— i6)) = Sy(£ = —z(t))\[k]) determining "principal" parts in the 
evolution operators Un(t) = exp(ia(t)Y )SY(£(t)) for initial states |0(O])) = |[/j])[13]. 
Similarly, this approach yields nonlinear quasiclassical Bloch-type equations 

y = ^x V C, y = (2/1,2/2,2/0), C = ±y 2 +y 2 1 +y 2 2 , V = (d/dy u d/dy 2 ,d/dy ) (3.17) 

which are equivalent in the mean-field approximation (3.14) to those obtained from 
Eqs. (3.7) (yi = z\{Yj}\ [Zj]; z), A x B is the vector product symbol). 

So, Eqs. (3.11)-(3.17) yield 5X(2)-cluster quasiclassical solutions of spectral and 
evolution problems which take into account quantum correlations of interacting sub- 
systems though they do not describe quantum dynamics exactly unlike the case of 
s/(2)-linear Hamiltonians [13]. 



4 . POLYNOMIAL OSCILLATOR LIE ALGEBRAS in ANALYSIS of COM- 
POSITE FIELD MODELS 



Another natural area of applications of algebras g p a is an algebraic analysis [5,18] 
of composite field models with internal (gauge) symmetries [2] which generalizes basic 
ideas of the paraquantization [3]. The simplest example of such analysis (but without 
introducing algebras g p( j) was given in [18] by using models (3.2) with m = = u 
to describe resonance states in particle physics; later it was generalized on multimode 
cases and applied to study multiphoton processes in quantum optics [5,10]. 

It was shown that operators V + = (af) n describe C ra -invariant n-particle kinematic 
clusters which display unusual (para) statistics and correspond to generalized asymp- 
totically free fields realized on the Fock spaces L F . Specifically, l^-clusters satisfy 
non-canonical CRs (3.1) and multi-linear relations [5]: 

adp' 1 V + = 0, ad v V + = [V, V% V = (ai) n = (V + ) + (4.1) 

generalizing (for n > 3) trilinear parastatistical Green's relations [3]. Furthermore, the 
subspaces L([l = k/ti]) in (2.2) describe coherent mixtures of constant numbers k of 
uncoupled particles and of varying in time numbers iVy of V-clusters. However, 
operators Ny have not standard (for (papra) fields) bilinear in V,V + forms (cf. [2,3]) 
but they can be expressed as nonlinear functions in the bilineals V + V, VV + [5]: 

N v = (E 11 -C(R ))/n, E n = ata 1 = nV = <p(V + V), C([l }) = nl = k (4.2) 

as it follows from Eqs. (3.1), (3.5). Therefore, at best the quantities V + can be 
set in correspondence only to parafield quanta [3] (when n — 2) rather than to any 
asymptotically free particles [5]. Nevertheless, one can construct operators W + = 
W + ({Vi}),W = (W + ) + obeying canonical CRs [W, W + ] = 1 and having the standard 



number operators N w = W + W(= N v ). Specifically, in [5,18,12] two equivalent forms 
were found for W + : 



Vo ~ Ro + 1 
(£ii + n)W 



,1/2 



w + = v + c r {v + y\vy = v 



+ 



(4.3) 



r>0 



where the second form is, in fact, a modification of the (inverted) map (3.10) [12]. 

The analysis above was generalized in [5] by means of: 1) extending on the case of 
"to" modes using C n -invariant Hamiltonians 



# = c + + £ [gn.vl+giy h j,v+.. = <•••<; (4-4) 



2) involving both boson (a^) and fermion (6j) variables; 3) considering Hamiltonians 
with non-Abelian invariance groups Gi = SU (n) (obtained from (4.4) by the substitu- 



where polynomial oscillator (super) algebras osc v (m;n) and osc x (m; 1") are extensions 
of the unitary algebras u(m) = Span{Eij = af aj} and Span{Eij = (aj 1- • aj)} by their 
symmetric (V^,V..) and skew-symmetric (X+,X_) tensor operators. 

The operators X+,X_ and V+, V„, satisfy non-canonical CRs with right sides de- 
pending on Eij (and on the SU(n) Casimir operators for osc x (m; 1™)) and obey (due 
to the invariant theory) certain extra "bootstrap" relations (of the type V1...1V2...2 = 
V21...1V12...2) [5] which are similar to those occuring in quantum field theories with con- 
straints [2,3]. All this entails unusual statistical and other properties of Gj-invariant 
clusters associated with X+,V+ and complicates extensions of the 1-mode analysis 
above. Specifically, the task of obtaining m-mode generalizations of the map (4.3), 



W+ = £ {VtJXt.JfUSiEij}), [W a , W b + ] = 5 ab , W a = (Wt) + (4.5) 



(with coefficients /"(•••) determined from finite-difference equations due to [W a , W^} = 
S a b) is, in general, sufficiently difficult because of relations abovementioned between 
VyX-clusters. 



In conclusion we briefly point out some of prospects of developing results obtained. 
For instance, formal constructions of Sections 2, 4 may be generalized by involving into 
consideration q-deformed oscillators and other invariance groups Gi. It is also of interest 
to examine infinite-dimensional algebras g p d related to g p d along lines of standard studies 
[7,8] and to construct non-Fock realizations of the g pa - IRs (cf. [11]) extending arbitrary 
/i-modules with the help of coset generators. 

It is of very importance to develop exact methods of Section 3 since they outline 
ways of generalizing standard group-theoretical technigues for solving both spectral and 
evolution tasks with dynamic symmetry algebras g p d- Specifically, one may use tech- 
niques of g-deformed algebras and g-special functions [6] (due to interrelations between 
sl pa -(2) and g-deformed algebras) as well as pseudogroup and braided geometry concepts 
[9,10] (due to Eqs. (3. 8), (3. 9)) for determining new classes of special functions related 
to g p d- At the same time quasiclassical approximations obtained may be considered as 
specific asymptotics of exact solutions that opens a possibility to use the techniques 
of asymptotic expansions [19] for finding latters. On the other hand, solutions of the 



m 



i=l 




5. CONCLUSION 



nonlinear Eqs. (3.7) enable to determine operator analogs of Abelian functions which 
are, probably, related to quantization problems on algebraic varieties [12]. 

Results of Section 4 provide an effective analysis of composite models with internal 
(^-symmetries at algebraic and quasi-particle levels, when obtaining explicit expressions 
for /...(. . .) in Eqs. (4.5) and examining the limit "m — > oo". Furthermore, involving 
into consideration spatiotemporal symmetries, one can construct appropriate "physical" 
composite fields and relevant Hamiltonians (Lagrangians) for them in terms of " quanta" 
W a (cf. [2,3,5,10]). 

The author thanks A. Odzijewicz for interest in the work, useful discussions and 
remarks and V.N. Tolstoy for stimulating discussions. A partial support of the Russian 
Foundation for Basic Research (under grant No 96-02 18746a) is acknowledged. 
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